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Starting from the Chiral Perturbation Theory Lagrangian, but keeping differ- 
ent masses for the charged and neutral mesons (m u 7^ mj), and using a pre- 
viously developed non-perturbative unitary scheme that generates the lightest 
meson-meson resonances, we construct KK — > KK and KK — > tt + tt~ in the 
vector channel. This allows us to obtain the kaon-loop contribution to the (f> — p 
mixing and study the <j> — > it + it~ decay. The dominant contribution to this 
decay comes from the 0^7^ 7r + 7r~ process. However, there can be large 
interferences with the subdominant contributions coming from <f> — p and <f>~Lo 
mixing, or of these two contributions among themselves. As a consequence, 
a reliable measurement of <j> —* n + ir~ decay could be used to differentiate 
between some cj> — u> mixing scenarios proposed in the literature. 

PACS numbers: 13.25.Jx, 14.40.Cs, 12.39.Fe 



I. INTRODUCTION 

The <j) decay into 7r + 7r~ is an example 
of isospin violation, since the <f> has isospin 
1 = and spin J = 1, and it would not 
couple to 7r + 7r~ in the isospin limit, which 
requires I + J = even (the decay into 7r°7r° 
is forbidden in any case because the parti- 
cles are identical). In addition, it violates 
the OZI rule jl| and hence it is subleading 
in the large N c M expansion. The experi- 
mental situation on this decay is rather con- 
fusing. There are two old results: BR = 
(1.94 + 1.03 - 0.81) x 10~ 4 from J|, and 
BR = (0.63 + 0.37 - 0.28) x 10~ 4 {j, with 
very different central values but whose er- 
rors are so big to make them compatible. 
Very recently, two new, more precise, but 
conflicting results have been reported from 
the two experiments at the VEPP-2M in 
Novosibirsk: the CMD-2 Collaboration re- 
ports BR = (2.20 ± 0.25 ± 0.20) x 10~ 4 § 
whereas the SND Collaboration || obtains 
BR = (0.71 ± 0.11 ± 0.09) x 10~ 4 . 

On the theoretical side, the common 
ground is based on the 4> — p mixing 
to account for the strong part of the decay. 
In addition, in ref. || it has been pointed out 



that the two-step <j> — lo — p transition^ can 
give a relevant contribution and that other 
non-resonant processes, as a possible bare 
4>pn coupling, have to be considered in de- 
tail. It is remarkable, in contrast with the 
OZI allowed u> — *■ tc + it~ decay, that the elec- 
tromagnetic </i>7r + 7r~ coupling via photon ex- 
change (f> — j — p — tt + tt~ provides the right 
order of magnitude [^J|] . 

Within Chiral Perturbation 

Theory (ChPT) |o]Q, isospin breaking has 
recently gained interest, since it is possible 
to take systematically into account the cor- 
rections due to the different u and d quark 
masses and to electromagnetic effects. Ex- 
amples of such calculations are 7T7t scattering 
|l2j , some ttN amplitudes and the nucleon 
self-energy p3| ; NN scattering and the 
pionium atom 

Unfortunately, isospin violation in <j) ~ * 
7r + 7r~ lies far away from the ChPT appli- 
cability range, since it involves the propa- 



* As a matter of fact, this two-step process, just 
gives a contribution to the (j> — p mixing. We will 
consider such resonant processs as the one that 
provides, by resonance saturation, the comple- 
mentary local terms to the kaon loop contribu- 
tions to <j> — p mixing that we will calculate later 
on. 
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gation of the pair of mesons around 1 GeV. 
Nevertheless, new nonperturbative schemes 
imposing unitarity and still using the ChPT 
Lagrangian have emerged enlarging the con- 
vergence of the chiral expansion JTq-p^l , for 
a review see ref. Here we shall fol- 

low the work |TJj] , since it provides the most 
comprehensive study of the different meson- 
meson scattering channels, including reso- 
nances up to 1.2 GeV. In particular, this 
method yields a resonance in the / = 0, J = 
1 channel, which is related to the </> and thus 
will allow us to obtain an important con- 
tribution to 4> ~~ * 7T + 7r~ due to the charged 
and neutral meson mass difference. We shall 
also consider electromagnetic contributions 
at tree level as well as the contribution due 
to the <f> — uj mixing. These three contri- 
butions can have different kinds of cancella- 
tions among themselves, depending on the 
<f> — uj mixing scenario. 

Some other theoretical uncertainties in 
our approach are unavoidable since the re- 
sults are rather sensitive to the Li coeffi- 
cients of the 0(p 4 ) ChPT Lagrangian and 
to the value of F v , which measures the cou- 
pling of a vector resonance with a photon. 
We will not calculate the electromagnetic 
loop corrections since the present ignorance 
of higher order counterterms makes their 
calculation unfeasible. However, from refs. 
|7|, p0|j2lJ| one expects the meson-photon in- 
termediate states to yield a contribution of, 
at most, 25% of that of kaon loops. 



II. TREE LEVEL CONTRIBUTIONS 

A. The vector meson chiral Lagrangian 

In order to calculate the contribution of 
an intermediate photon to cf> — > it + tt~, we 
will use the vector meson chiral effective La- 
grangian presented in [ p2[ , which is written 
in terms of the SU(3) pseudoscalar meson 
matrix (j) and the antisymmetric vector ten- 
sor field VJj„ defined as 



The latter is normalized such that 
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-[P^(R) - P^(R)}, (2) 



with Mr, P and e^(-R) the mass, momen- 
tum, and polarization vector of the vector 
field R. Following |22| , let us then consider 
the Lagrangian 

^[v(i--)\ = ^<wr> + l -^{v^ U "), 

(3) 

where "( )" indicates the SU(3) trace and 



= iv)D^JJu + = uj t 




DJJ = d tl U-ie[Q,U]A tl , (4) 
with Q the quark charge matrix 

Q = ±diag(2,-l,-l), (5) 

and the electromagnetic field. As usual, 
/ is the pion decay constant in the chiral 
limit (we take / ~ /„. = 92.4 MeV) and the 
fft u and F^ v tensors are defined as 



F» v = eQ(& l A v -d v AV), 



(6) 



In order to introduce the <fi an d w states, 
we extend SU(3) to U(3) and substitute 



(7) 



where -Z3 is the diagonal 3x3 matrix and 
u>i is the lightest singlet vector resonance. 
Hence, by imposing ideal mixing between u>i 
and us 
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Unless otherwise stated, in the following we 
will refer to these states simply as uj and </>, 
although it should be kept in mind that we 
are referring to their respective ideal states. 
Finally, the <j) and uj can be introduced into 
the chiral notation by replacing in eq. (|^) the 
tensor by 
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The convention of signs of eq.({|) agrees 
with a more standard one if we take e neg- 
ative in all the Lagrangians, as we shall do 
in what follows. The vertex function 4> — > 7, 
resulting from eq.(0), is 



andeq.ljjj), substituting g^x+K- by 5 07r+7r - 
and px + by , provides the tree level elec- 
tromagnetic contribution to the (f> — > 7r + 7r~ 
decay width. With a value of F v = 154 MeV 
from the p — * e + e~ decay |2^] this contribu- 
tion alone would yield BR(<p — >■ 7r + 7r~) = 
1.7 x 10 -4 , a value compatible with the ex- 
periment of ref. Q, within errors. 



B. Comparison with ui 
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and to the same order than eq.(|To|) the La- 
grangian giving the coupling of the photon 
to the pions is 
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With these ingredients we can write the 
contribution of the Feynman diagram of 
fig.l, which is given by 
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n<f)(p«+-P«-)»F{Ml), (12) 



It may seem surprising that 9^j+ n - al- 
ready provides the correct order of magni- 
tude of the <p — *■ 7r + 7r~ decay, since, in con- 
trast, it is well known 25 2(| that the tree 
level photon contribution uj — > 7 — > p — ► 
7r + 7r~, represents a negligible amount of the 

T(UJ — > 7T + 7T~). 

However, the case of the w is radically dif- 
ferent from ours and can be well understood 
from p — uj mixing. We will now calculate 
this effect making use of an effective chiral 
Lagrangian and large N c arguments |27| . In- 
deed, from this reference, the p — uj mixing 
can be represented as 



(16) 



where F{q 2 ) is the pion electromagnetic 
form factor, which at the <f> mass is given 
by F{Ml) = -1.56 + i 0.66 §|. 

This can be compared with the coupling 
of the 4> to K + K~, or K° K°, which can 
be obtained from the Gy term in eq. (|3|) and 
reads 

i£<j>K+K- -> -ig^K+K- e M (0) (PK+ ~PK-)ti 7 
sG v 

9 * K+K ~ ~ V2fM~i 
The 4> — > K + K~ width is then given by 



r — Pk+ rfi 

L4>K+K- ~ g7rM 2 9<PK+K- 
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which, using its experimental value p4) , pro- 
vides Gy — 54.3 MeV (to compare with 
Gy = 53 MeV. from the study of the pion 
EM radius @,g2|). 

By analogy to eq. (|l3"|) , eq.(p^) provides a 
(j>TT + n~ coupling 
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Note that My is the mass of the vector octet 
in the chiral limit M v &_M P [0 . We have 



made use in eq. (|l6| ) of eq. (|2|), thus 
turning to the usual vector notation. At low- 
est order in ChPT [jllj the first two terms 
in eq. (|l7]) arise from the quark mass dif- 
ference, and the third one is of electromag- 
netic origin from the exchange of a photon 
between the p and the uj. It is straightfor- 
ward to see that the electromagnetic contri- 
bution only amounts to a 14% of that due 
to quark mass differences. 

Contrary to the <j> case, the p — uj mixing is 
OZI allowed and leading in large N c , as can 
be seen from eqs. (|l6|) and flT?|). In fact, 
this term is of the same order than the free 
Lagrangian, both in the l/N c and in chiral 
countings (this is more clearly seen in the 
tensor notation). 

In addition, there is a kaon loop contribu- 
tion, fig. 2, which, from ChPT, is expected 
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to be of the same order of magnitude than 
the electromagnetic contribution. Evaluat- 
ing the diagram of fig. 2 one has: 



c 2 s 2 
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[L(s,m K +) - L(s, 
(18) 



where, once again, we have used eq. (||) 
to present our results in the vector nota- 
tion. The L(s, m) loop function, in the usual 
ChPT ITS - 1 scheme, is: 



the u> — ► ir + ir~ decay comparatively much 
larger. Note that such contribution is miss- 
ing in the OZI violating, large N c sublead- 
ing, 4> — P mixing. The fact that to — > 7r + 7r~ 

is much larger than the d> — ► 7r + 7r~ dominant 
m K o)X. . to . / . 

contribution is very relevant since, through 

the 4> — lu mixing, it provides an additional 

mechanism that has to be taken into account 

in the complete calculation of <j) — > 7r + 7r~, 

and that we analyze next. 



C. The "two step" 
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As a matter of fact, the physical </> and 
lu states are not the ideal ones defined in 
eq.(ph, but instead 
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where p, is the dimensional regularization 
scale. In order to estimate the eq. ( |I§| ) 
contribution at s = M 2 , we use the natu- 
ral value p = A^ n prp « M p . The results 
depend on the regularization scale but they 
provide a good estimate of the order of mag- 
nitude, as we shall see later on, when we will 
reevaluate this contribution within the chi- 
ral unitary approach. 

At this point we are ready to compare all 
contributions: 



10 ~ k)—™' _ 8y<f) 1 ' 

4> ~ 5 V io ldeal + 4> ldeal 

In the literature there is a general agreement 
on \Sv \ — 0.05, but, apart from conventions, 
not on its sign ||. Its contribution to the 
4>Tr + 7r~ effective coupling is obtained from 
fig. 3 as follows: 
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Quark mass differences from eq.(^) =— 5221. SWePva^e already obtained 9 pu , although 
EM contribution from eq.(|l7|) = 725.h<Me^i 2 has to be evaluated at \fs — M<p. 

Kaon loops from eq.([T^) = —13(SffiIb , \tne dominant contribution comes from 

the quark mass differences. The kaon 



Hence, the p — to mixing is dominated 
by the OZI allowed strong contribution due 
to quark mass differences, which is leading 
both in the large N c and chiral countings. 
In addition, the kaon loops are smaller than 
the electromagnetic contribution although 
with a large destructive interference between 
them (for Gy = 65 MeV, which is the value 
needed to reproduce T(p — > 7r + 7r~) from eq. 
(||), the estimate of the kaon loop contribu- 
tion would be -190 MeV 2 ). We will find 
again this large destructive interference be- 
tween the kaon loops and the electromag- 
netic contribution when considering the </> 
resonance. 

In summary, the fact that the purely elec- 
tromagnetic contribution already provides 
a reasonable order of magnitude for the 
<f> — > 7r + 7r - decay, is due to the absence of 
the OZI allowed contribution, which makes 



loop contribution cannot be calculated us- 
ing eq.(|i~8|), since that formula is not unitary. 
We will see later, how this number can be 
obtained from the chiral unitary approach, 
and again it is of the order of 200 MeV 2 , and 
therefore numerically irrelevant for the fol- 
lowing discussion. 

The new parameter can be obtained 
from the literature. Nevertheless, its imag- 
inary part can be obtained from unitarity. 
The most relevant intermediate states are 
K K and three pions. In the first case the 
couplings to <f> and to are completely deter- 
mined by the vector resonance Lagrangian. 
However, the imaginary part contribution 
of three pion intermediate states has some 
model dependence ||, mostly through the 
g^pn coupling. 

We consider now two different scenarios 
for the <f>— io mixing which illustrates to some 
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extent the uncertainties that are found in 
the literature with respect to this issue: 

• "Weak mixing" scenario ||, where 
Re e«ju, = and = 0.78 GeV -1 . 

• "Strong mixing" scenario ||, where 
ReO^ = 20000 to 29000 MeV 2 and 
94> P tt = 0. 

which will therefore appear as different cases 
in our final result. 

Up to now we have just concentrated on 
the tree level diagrams of the </>7r + 7r~ decay. 
There are, however, important contributions 
from kaon loops that we will analyze in the 
next sections, whose calculation is the main 
novelty of this work. 

III. DIRECT KAON LOOP 
CONTRIBUTION TO p MIXING 

A. Introduction 

The pure strong interaction chiral La- 
grangian gives a contribution to the <f> — > 
7r + 7r~ decay if the charged and neutral me- 
son masses are different, otherwise it would 
be forbidden. 

For instance, from eqs.(^) and ([)]) there 
is no direct 07r + 7r~ coupling. However, we 
can generate a non vanishing <f> — > tt + tt~ 
transition when keeping different masses for 
the charged and neutral kaons in the loops 
of fig. 4, which do not violate the OZI rule, 
although they are subleading in large iV c . 
In fact, these diagrams are expected to give 
the main strong interaction contribution to 
4> — > 7r + 7r~ due to intermediate states. For 
instance, the 4> couples much more strongly 
to K K than to 3tt, as it is clear from the fact 
that T(> -> 3n)/T(KK) w 1/5, although 
three pions are kincmatically much more fa- 
vored than two kaons Q 

Note that in the evaluation of the dia- 
grams of fig. 4 the KK — > ir + ir~ amplitude 
can receive important contributions from 
the u> or p exchange. In the first case, the 



"We will address in subsection III.D the prob- 
lem raised in refs. | p7|4q ], relative to the con- 
tributions of more massive virtual intermediate 
states 



uj couples to the p once again, and there- 
fore is included in the tf> — lu — p mixing con- 
tributions. Thus, in the following we will 
concentrate on the evaluation of these kaon- 
loop contributions to the direct 4>— p mixing, 
that is, we will consider only the exchange 
of the p in the KK — > n + ir~ 1 = P-wave 
amplitudes appearing on fig. 4. 

An estimation of the imaginary part of 
this contribution to the diagrams in fig. 4 is 
straightforward using the vector meson chi- 
ral Lagrangian. The sum of the diagrams 
does not vanish due to the different masses 
of the charged and neutral kaons. The 
cf> — > 7T + 7r~ branching ratio that would be 
obtained taking into account just this con- 
tribution is already of the order of magni- 
tude of the experimental results, given the 
large uncertainty on the data. 

Yet, this estimate does not take into ac- 
count corrections to KK — > n + ir~ due to 
isospin violation. In addition, the real part 
of the loop remains ambiguous since it re- 
quires the knowledge of higher order con- 
tributions than those given by eq.(||), that 
is, counterterms to absorb loop divergences. 
Furthermore, even when we have such coun- 
terterms, the chiral expansion is only ex- 
pected to work at energies which are below 
the <fi mass. 

B. Resonances and the I AM 

We present here a method which deals si- 
multaneously with all these problems in or- 
der to extract the aforementioned kaon loop 
contributions. The method exploits the in- 
formation of ChPT up 0(p 4 ), by relying 
on the expansion of the T~ ^matrix. The 
technique starts from the 0(p 2 ) and 0{p A ) 
ChPT Lagrangian and uses the inverse am- 
plitude method (IAM) in coupled channels. 
Unitarity provides for free the imaginary 
part of T , and then a chiral expansion 
is done for ReT -1 , which, in the present 
case, has a larger radius of convergence than 
T itself. This approach has been applied 
in the isospin limit with remarkable results: 
with just one channel |l6| it nicely describes 
the cr, p and K* regions, amongst others, 
in 7r + 7r~ and 7rA" scattering. When gener- 
alized to coupled channels |l7],[l8| it also de- 
scribes meson-meson scattering with all the 
associated resonances up to about 1.2 GeV. 
A more general approach is used in [ p8| by 
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means of the N/D method, in order to in- 
clude the exchange of some preexisting reso- 
nances explicitly, which are then responsible 
for the values of the fourth order chiral pa- 
rameters. 

The T amplitude is defined in terms of the 
partial waves as 

T = ^(2J+l)T / (s)P J (cos0), (21) 
J 

In what follows we will refer to Tj simply as 
T. Within the coupled channel formalism, 
the IAM partial wave amplitude is given by 
the matrix equation 



T = T 2 [T 2 -T i ]- 1 T 2l 



(22) 



where T 2 and T 4 are 0{p 2 ) and 0(p 4 ) ChPT 
partial waves, respectively. In principle, T4 
would require a full one-loop calculation, but 
it was shown in [ p"7| that, at the phenomeno- 
logical level, it can be well approximated by 



Re T 4 ~ T 4 + T 2 ReG T 2 



(23) 



where T± is the tree level polynomial 
contribution coming from the £4 chiral 
Lagrangian and G is a diagonal matrix 
diag(gi, g 2l 33), where gi is the loop function 
of the intermediate two meson propagators, 
which we give in the appendix. In JlTJ the 
loop integrals are regularized by means of a 
momentum cut-off, q m axt m the loop three- 
momentum. The relation between this cut- 
off and the dimensional regularization scale 
/i, normally used in ChPT, is also given in 
that paper. 

We have also taken advantage to correct a 
small error detected in |L7| in the K + K~ — > 
K°K° amplitude, whose complete expres- 
sion in the isospin limit is given in the ap- 
pendix. We have also reconducted a fit to 
the data including those on (Sqq—Sh), which 
are well determined from |29|. The fit of 



the phase shifts and inelasticities is carried 
out here in the isospin limit, as done in p7[ . 
There are several sets of Li coefficients which 
give rise to equally acceptable fits. 

As it can be seen in fig. 5, there are sev- 
eral plots for which there are incompatible 
sets of data. This is particularly evident for 
the <5oo data both in mr — > mr and KK, in 
the inelasticity 7700, and in the <5oi/2 phase 
shifts. As a consequence, although we have 
performed a x 2 fit of the data using MINUIT 
pOf , the resulting x 2 P ei ' degree of freedom is 



not really very meaningful, since the Li val- 
ues depend on the estimate of the systematic 
error of each experiment, which is not given 
in many original references. In addition, due 
to the fact that we have eight parameters, 
there are several x 2 minima, which yield 
very similar values of x 2 f° r rather different 
values of some chiral parameters. Which one 
is the real minimum depends on how we add 
the systematics. For that reason we have 
preferred to give several sets of coefficients, 
which yield x 2 ''/d.o.f < 2 when assuming a 
3% systematic error added in quadrature to 
the statistical error quoted by each experi- 
ment. 

We write in table 1 the different sets of 
chiral parameters and we show their cor- 
responding results for the phase shifts and 
inelasticities in fig. 5. We can see that the 
small differences in the results appear basi- 
cally only in the a (980) and k(900) reso- 
nance regions, where data have also larger 
errors or are very scarce. 

Although the tadpoles and loop terms in 
the crossed channels were neglected and re- 
absorbed into Li redefinitions |[7) when we 
use eq.(|2^), these coefficients are still close 
to those of standard ChPT (see Table I). 
Consequently, it seems that this simplify- 
ing approximation has a small effect in the 
relevant energy region, not spoiling also the 
standard low energy ChPT results. 

One of the side consequences of the ap- 
proach was the generation of a resonance 
around 1 GeV in the 7 = and J = 1 chan- 
nel, which only couples to KK. Actually, it 
has a zero width, since its mass is below the 
KK threshold. One is tempted to associate 
this state to the <f> meson, however, we can 
only relate it with the octet part wg, which, 
by mixing with a singlet generates the </) and 
the lu. This can be easily understood since 
the singlet in this channel, u>i , which is sym- 
metric in the SU(3) representation, does not 
couple to two mesons because their spatial 
wave function is antisymmetric. Since only 
two meson states were considered in ]T?j^q , 
0J1 does not appear in the IAM, and the res- 
onant state found in that channel can only 
be related to ug. However, we will see next 
that we can still exploit the properties of the 
wg pole in order to study the decays of the 



resonance. 
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C. Extracting the <f>n tt coupling from 
the IAM 

Let us then turn to the case of interest 
for this work: the evaluation of the J = 1 
KK — > KK and KK — > it + ii~ amplitudes 
around the mass of the u)%. Now we are 
breaking isospin explicitly by keeping differ- 
ent the charged and neutral meson masses, 
while keeping the Li obtained from the pre- 
vious fits to meson-meson scattering in the 
isospin limit. In addition, we are dealing 
with three two-meson states: K + K~ , K°K° 
and tt + it~, that we will call 1, 2 and 3, re- 
spectively. The amplitude is a 3 x 3 ma- 
trix whose elements we will denote as Tjj 
(for instance, T13 stands for the J = 1 
K+R- -> n+n- amplitude). The T 2 and 
Tf amplitudes used in the present work and 
calculated in the isospin breaking case, are 
collected in the appendix. 

Once the amplitudes are unitarized with 
the IAM, one observes the presence of two 
poles, one corresponding to the p(770) and 
the other one to the w% resonance. It is inter- 
esting to note that the w$ pole appears with 
a mass around 910 MeV, very close to the 
value 930 MeV predicted by the quadratic or 
linear SU(3) mass formulae jhj| for the Wg 
mass. In the following, we will denote by fig 
the resonance pole that we have obtained in 
our approach corresponding to the wg reso- 
nance. The motivation for this change of no- 
tation is the lack of the 37r state in our model 
since this contribution can be particularly 
relevant in order to study certain properties 
of the u>g resonance. For instance, the 3-7T 
couplings of the u>g and w\ according to eq. 
(g) add for the uj giving rise to the w — > 3-7T 
coupling and almost cancel each other in the 
case of the <j>, |fl^->3jr| << ISw— 3tt|- 

In order to evaluate the kaon loop con- 
tribution to the (j)TT + ir~ coupling via direct 
4> — p mixing, we first study the fig7r + 7r~ 
coupling. We thus evaluate the K + K~ — > 
K+R- amplitude (T u ) and the K + K~ -> 
7r + 7r~ amplitude (T13) near the pole of the 
fig resonance. Close to the i7§ pole the am- 
plitudes obtained numerically are then dom- 
inated by the exchange of this resonance, 
represented diagrammatically in fig 6. 

By considering couplings like those in 
eq.((l|) for the fig to K+K~ and tt+tt", 
these two amplitudes, once projected in the 
,/ = 1 channel, and close to the fig pole, are 
given by 
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s lvl n 8 6 

where pi is the modulus of the center of mass 
three-momentum of the i particle. The dia- 
gram of fig. 6b can be interpreted as provid- 
ing an effective strong % „.+„.- coupling. 

By looking at the residues of the Tu and 
T13 amplitudes in the pole we can get 

9q s k+k- 9n s K+K~ an( i 9q 8 k+k- gn^+n- ■ 
Thus, defining 



= lim 



we obta 



9q. 8 k+k- 



Ql3 

Qu 



(25) 



(26) 



In eq.(|25|) the Tn,Ti 3 amplitudes have a 
large p exchange background, which can be 
eliminated using the residue of the fig pole 
obtained via eq. (|25|) . Yet, numerically this 
background can be eliminated to a large ex- 
tend by using the isospin zero combination 
-(K+K- + K°K°)/V2 in the initial state. 
Hence, the <7n S 7r+ir- is more efficiently eval- 
uated by means of the combination 



9n 8 K+K- 



Q13 + Q 



2:i 



(27) 



We have checked numerically that the 
9si 8 k+k- and the g n8 K°R° couplings have 
the same value in our approach. Since we 
are interested in (f>, we still have to make 
the connection between the fig7r + 7r~ and 
(/>7r + 7r~ couplings. Indeed, we have explic- 
itly checked that, when removing the rescat- 
tering resummation implicit in the IAM (by 
setting G — 0, see eq. (f23|)), the ratio in eq. 
fl27| ) becomes between one and two orders 
of magnitude smaller. Even more, this dras- 
tic reduction in the fig7r + 7r~ coupling is also 
obtained when making G = diag(0,0, gs), 
that is, when only removing the kaon loops. 
Therefore the fig decays to ir + ir~ mainly 
through the mechanism shown in fig. 4 (re- 
placing the (f> by the fig). This observa- 
tion allows us to find the kaon loop contri- 
bution to <p — > 7r + 7r~ that we are looking 
for, through the same mechanisms of the fig, 
since the only difference will be the initial 
tt$KK and tpKK couplings, which can be 
canceled taking the following ratio 
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9(j>K+K- 



Therefore, from eq. fl2q), one has 



ffn 8 7T+7T- 

9n s K+K- 



(28) 



( S ) / V Ql3 + Q23 , x 

' 1 ~ n _Ln 9<t>K+K-{s), ( 29 ) 



with g^K+K- given in eq. Jl3|). Here we 
are neglecting the mass difference between 
the fig and the 4> resonance, which is around 
100 MeV. In any case one has to take into ac- 
count that: 1) The important p exchange ef- 
fect is also canceled in the ratios. 2) We have 
removed in eq. ( p5[ ) the three-momenta fac- 
tors. As a result, the remaining differences 
coming from the mass difference should be 
rather tiny. 

Finally, by adding the above contribution 
with that of eq. (|l5|) , we find 



(7) (4>u>) (s) 



(30) 



which allows us to obtain the <f> — ► tt + tt 
decay width as we did before only for the 
fldX+ir- coupling. In order to determine the 
sign of the interference in eq. ( |30"| ) it is impor- 
tant to know the sign of Fy Gy (see eqs.([l3|) 
and (p^)). We have taken Fy Gy > since 
the Lg chiral parameter, whose main reso- 
nance contribution is given by Fy Gy /2M^ 
p2[, is positive and large. 



D. The OZI rule violation 

The direct coupling gs-K+n- violates the 
OZI rule. This is clearly seen in a quark 
picture when considering the <f> as a pure 
ss state. From the QCD Lagrangian one 
can see the OZI rule as a prediction of the 
1/N C expansion, with N c the number of col- 
ors. While the couplings of the decays which 

1 /2 

do not violate the OZI rule are 0(1/N C ) 
those that violate the OZI rule are sup- 
pressed by an extra 1/N C . In addition, me- 
son loops are suppressed by at least one 
power of l/iV c M. As a consequence, the 
g£+ n - coupling given in eq. @, which 
is due to kaon loops, as discussed above, 
is C(l/iV 3 / 2 ). Note, in contrast, that the 
94>k+k- coupling, from eq. (p"3|), is order 

l/Nc /2 , since / and Gy are 0{Nl /2 ) and 
Ms is order 1. 



However, in quark model calculations p7[ 
the large N c suppression of two intermediate 
meson states is considered insufficient in or- 
der to explain the experimental success of 
the OZI rule. The point is that in these 
models the real parts of the two meson loop 
contributions to OZI violating processes, al- 
though large N c subleading, are found to be 
much larger than they should be in order 
to explain the experimental success of the 
OZI rule. The solution advocated by the 
authors is that a cancellation among a very 
large number of intermediate states seems to 
operate. This is illustrated via the example 
of uj — p mixing in |47| . 

Nevertheless, one should notice that the 
real part of the two-meson loop is diver- 
gent, and the remnant finite part depends 
upon the regularization and renormalization 
schemes, apart, of course, from the details 
of the dynamical model. In refs. Jl7|,fl8) this 
regularization is done including several cut- 
offs within a quark flux tube model, having 
an explicit scale dependence. In contrast, 
we have just included kaons and pions as 
intermediate states and we have renormal- 
ized such contributions making use of a cut- 
off w ^ChPT' Still, the physical quanti- 
ties we calculate are scale independent and 
well defined, since any change in the cut- 
off would be reabsorbed by a change in the 
Li ChPT counterterms. Note that, since 
we are making use of an effective field the- 
ory formalism, the chiral Lagrangian coun- 
terterms should take into account any other 
contribution from more massive intermedi- 
ate states. In our approach we use ChPT up 
to Ofp 4 ) and generate higher orders through 
eq. (|22|). In this way, any other contribu- 
tion coming from heavier virtual interme- 
diate states is reabsorbed in the final val- 
ues of the Li counterterms given in table 1. 
At this point, our previous statement about 
the fact that our result for the q\ • , _ cou- 
pling is due to kaon loops is meaningful only 
because we have taken a natural value for 
the cut-off. For such value, the contribu- 
tion from graphs without kaon or pion loops, 
which come just from the Li counterterms, 
is between one and two orders of magnitude 
smaller than that of kaon loops. Comparing 
our work with that of refs. p7 48 1, we cannot 
tell exactly the size of each separate contri- 
bution due to each state more massive than 
the kaons. If each one of these contributions 
was large as it happens in refs. |47|]l8j| , then, 
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we would also be finding a cancellation. 

In order to obtain further support for our 
arguments about the kaon loop size, it is 
instructive to revisit, within the IAM for- 
malism, the kaon loop contribution to uj — > 
7r + 7r~ that we estimated in sect. II. B. Note 
that the value obtained for the u> — p mix- 
ing from kaon loops in sect. II. B was depen- 
dent on the regularization scale. In con- 
trast, in the IAM this dependence is can- 
celed with that of the chiral parameters 
Li. In addition, the IAM respects unitar- 
ity and accounts for isospin breaking not 
only in the loops (through different masses 
of the charged and neutral kaons), but also 
in the 4>K + K~ and <f>K°K a couplings and 
the KK — ^ tt + tt~ amplitudes. 

In order to reinterpret our results for the 
f2g7r + 7r - coupling in terms of an Q% — p mix- 
ing and compare with sect. II. A., we write 
(see fig. 7) 



5n 8 



®n 8P g p ^ 



1 



ML 



(31) 



with g pn+n - = -G v s/(f 2 M 2 p ) from eq.(§). 

This gives us ©si 8P , from where, using eq.(||) 
and the fact that the uji does not couple to 
KK at the leading chiral order, we obtain 



1 9n 8 7T+7T- 

a/3 5p7r+7r _ 



[Ml 



(32) 



Taking now the value for <?n 8 7r+7r- obtained 
in the IAM from eq.(|27]), with gn s K+K- = 
— ^/S/lg^K+K- from eq.(||), we arrive at a 
value of Q U p{M p ) = (-52 - i76)MeV 2 and 
e wp (M ) = (-299 - z81)MeV 2 . These re- 
sults corroborate the "order of magnitude" 
arguments given in sect. II. B., obtained us- 
ing the non-unitary eq.(|l8|), to show that 
the kaon loop contributions are very small 
relative to the dominant OZI allowed con- 
tribution. 

It is also interesting to remark that the 
cancellation between mesons loops in the 
model of ref. ^] does not operate for the 
scalar sector with vacuum quantum num- 
bers J PC = ++ as discussed in ref. Q. 
The failure of the large N c suppression in 
this sector, and its associated OZI rule vi- 
olation, is also discussed in more general 
terms in ref. |49j]. Although the scalar 



sector is very hard to discuss in terms 
of quark models, due to the large rescat- 
tering effects, it is equally well described 
as the vector channels in the framework 
of non-perturbative unitarity methods from 
the ChPT series 0|J||g£l), see also 
fig. 5. For instance, in refs. J50| , |5l[ the a, 
/o(980) and a (980) were dynamically gen- 
erated and their meson-meson and 77 de- 
cay modes were analyzed in very good agree- 
ment with experiment. Furthermore, in ref. 
p8t the spectrum in the scalar sector was 
discussed taking into account as well the 
large N c limit. In addition, the presence of 
a scalar nonet due to the meson-meson self 
interactions, which disappears in the limit 
N c — > co, was then established. On the 
other hand, it was also found that the light- 
est preexisting scalar nonet, with mass 0(1) 
in the iV c counting, should comprise a sin- 
glet around 1 GeV and an octet around 1.4 
GeV, in qualitative agreement with the ex- 
pectations of ref. E^l • The success of our ap- 
proach in the + ~^sector indicates that our 
techniques are powerful in the study of OZI 
violating processes. Note that we describe 
both vector and scalar channels without in- 
cluding any new ad-hoc elements. 



IV. RESULTS AND DISCUSSION 



m 2 p +iM p t p ] a . 



4>- 



Branching Ratio 



In this section we are going to present the 
resulting branching ratios for the <j) — > ir + n~ 
decay. To do that we will consider and dis- 
cuss the different sources contributing to the 
total g^Tr+Tr- coupling as given in eq.(|30|). 

We first consider the contribution <?V +7r _ 
introduced in sec. II. A. We take as a fi- 
nal value ^ [10.6 ± 0.4 - i (4.47 ± 
0.15)]10~ 3 where the uncertainty is mainly 
due to the value of Fy, which ranges be- 
tween Fy = 154 MeV, coming from the 
p — > e + e~ decay, and Fy — 165 MeV, com- 



ing from the 



e + e decay, when evalu- 



ating both them with eqs. (||) and (||). 

Concerning the kaon-loop contributions to 
the 4>—p mixing eq. (|29|) , after averaging over 
all the fits presented in table 1, we obtain: 



0) 



-[5.6 ±0.4- i (3.8 ±0.12)] 10" 



Let us note that the error is mainly due to 
the differences between the Li correspond- 



9 



ing to the different fits, since they are much 
larger than the errors given by MINUIT, 
which are certainly underestimated. Fur- 
thermore, we have checked that this error 
band spans the dispersion in the results due 
to the variations of the chiral parameters 
that could yield a reasonable fit. 

Although they were not present in eq. ( |3C| ) 
there are corrections coming from diagrams 
with photon loops which are expected to be 
of the same order of magnitude as the isospin 
breaking corrections from the different mass 
of charged and neutral mesons p2[-^5[ . We 
do not have means at present to evaluate 
these diagrams within the non-perturbative 
chiral scheme which we have followed. One 
would also need counterterms whose values 
are nowadays unknown. However, explicit 
calculations of the absorptive part of the 777 
intermediate channel in ref. |?] give a con- 
tribution of, at most, 1/4 of the kaon loops 
but with opposite sign. This 777 will be our 
largest source of uncertainty in the errors 
given for each one of the different <f> — uj sce- 
narios, that we discuss next. 

As we have already commented, the con- 
tribution from the two step (f> — io — p mech- 
anism, depends on the <$> — uj mixing. Our 
results are the following: 

• Strong scenario: we find g'l^+ n - = 
[4.4- z3.7]10- 3 or g^_ = [6.0- 
«5.6]10~ 3 , depending on whether we 
use ReG^ = 20000 or 29000 MeV 2 , 
respectively. Therefore, there is a 
large cancellation with the kaon loop 
contribution, and we obtain: 

BR ~ (1.7 ± 0.3) x 10~ 4 to (2.5 ± 0.3) 

where the uncertainty in the central 
values depends on whether we use 
Re 6^ = 20000 or 29000 MeV 2 , re- 
spectively. 

• Weak scenario: we get <?^f +7r - = 
[-0.73 - i0.61]10~ 3 , very small com- 
pared with both the electromagnetic 
and kaon-loop contributions. Thus, 
there is only a partial cancellation of 
the electromagnetic contribution with 
that of kaon loops, and we obtain 

BR ~ (0.38 ±0.12) x 10~ 4 . 



Apart from the contributions discussed so 
far, there is also the possibility of local terms 
giving rise to a direct p — 4> mixing. How- 
ever, one can argue that, by resonance satu- 
ration, the inclusion of the two-step process 
(j) — to — p can be enough to take care of such 
local terms by considering that they are re- 
summed on the u propagator. 

V. CONCLUSIONS 

In this work we have evaluated the kaon 
loop contribution to the <fi — > n + Tr~ decay 
via (j> — p mixing from the splitting of me- 
son masses, making use of the unitarized chi- 
ral amplitudes with strong isospin breaking. 
We have shown that although this strong 
contribution to the — > 7r + 7r _ decay gives 
rise to smaller branching ratios by itself than 
the tree level electromagnetic contributions, 
they can have a very large destructive in- 
terference with either the electromagnetic or 
the <j> — uj — p contributions. 

We have also estimated the error in our 
4> — > 7r + 7r~ branching ratio calculation com- 
ing from the uncertainties in Fy, the fitted 
0(p 4 ) ChPT counterterms, the photon-loop 
contributions, as well as the considered 4>— u> 
mixing scenarios. 

A complete calculation of the loops with 
photons is missing in the present work, al- 
though they have been estimated making 
use of the results of ref. jjj. Still, they are 
the main source of uncertainty within each 
<p — u> mixing scenario. 

Accepting this additional uncertainty, we 
find that the strong coupling scenario Q 

x lW^ s 

BR ~ (1.7 ± 0.3) x 10~ 4 to (2.5 ± 0.3) x 10~ 4 , 

in very good agreement with the experimen- 
tal results of ref. 0. In contrast, the Weak 
§ scenario yields 

BR ~ (0.38 ±0.12) x 10~ 4 , 

It seems to prefer a value somewhat lower 
than the experimental value provided by 
ref. || , although still reasonably compatible 
with it. 

Of course, a precise determination of the 
photon loops in the non-perturbative regime 
would be desirable to reduce the theoretical 
uncertainties. 
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Finally, we would like to remark that the 
solution of the experimental conflict in the 
<f> — > 7r + 7r~ will, eventually, help us to dis- 
card some of the tp — uj mixing scenarios pro- 
posed in the literature. 



In the above formulas, /„-, f K + — 1.22 / OT 
and f K o arc the decay constants of the 
charged pion, kaon and neutral kaon, respec- 
tively. In the approach we are following here 
of neglecting tadpoles one has, up to 0(p 4 ), 
that 
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APPENDIX A: AMPLITUDES 

In this appendix we give the expression for 
the J — I partial waves obtained from the 
ChPT Lagrangian, but setting m u ^ m^. 
The normalization of the T-matrix used here 
is the same as in ref. ||l7|| . Let us first de- 
fine the modulus of the CM momenta of the 
different particles as 



P*+ 





4 



,Pk+ 



m' , . 

4 A + 



Pro 



where is the charged pion mass. Then, 
once they are projected in P-wave, the tree 
level amplitudes from the 0(p 2 ) and 0(p A ) 
Lagrangian for K + K~ — > 7r + 7r~ scattering 
are 



T 2 (s,t,u) = - 



Pir+PK+ 



T[{s,t,u) 



4 



[L 3 s 



For K+K- 
T 2 {s,t 1 u) 



K° 



K+ 



K+K 
2 



we obtain 



Pk+, 



)• 



T[(s,t,u) 



4p 2 



K+ 



3/ 4 



[2(2Li -L 2 + L 3 )s - 4(2L 4 + L 5 )m 2 K+ ] 



K+ 



the K°K° -> K°K° amplitude is exactly 
the same, but changing m K + by m^-o and 



Ik+ by f K °- For tt+tt 
2 



7T+7T , we find 



T 2 (s,t,u) 
T[(s,t,u) 



3/» 

*pU 

3/» 



[(2ii - L 2 + L 3 ) a - (4L 4 + 2L 5 )ml 



We have left the K+K~ -> K°K° ampli- 
tude for the end, since we had an erratum 
in our previous paper [fl7|| . Thus, we first 
give the complete amplitude in the isospin 
limit, before projecting on the P-wave. It 
reads 



T 2 (s,t,u) 



2m 2 K 



T 4 (s,t,u) = -J- [(4Li + L 3 )(s - 2m 2 ,) 2 + 
s „ Jk 

|2 



*K°- 



(2L 2 + L 3 )(t 



2m 2 K 



k(L. 



2um 2 Y Lz, — &m 2 K (2m 2 K 



8m 4 
- s)U] 



2L 2 (u - 2m 2 K ) 2 
2L 6 ) 



The P-wave in the isospin breaking case is 
given by: 



T 2 {s,t,u) 



T[(s,t,u) 



Pk+Pk<> 
3 Jk+ Ik° ' 

^Pk+Pk" r L 



3S 



L 5 (m 2 K -i 



rn 



whereas for K°K° 
are given by 



- L 5 (m 
7r + 7r~ scattering they 



3 f K + f K o 

)] PPhffisfiy^ we give the loop function G 
diag(g 1 ,g 2 ,g 3 ), where g t is 



771 



K 



»)] 



9i(s) 



(4tt) 2 



filog 



0~i Qi 



1 



o~i Qi — l 



2 log 



-{l + Qi) 



T 2 (s, t, u) 
T[(s,t,u) 



3 /if o f n 
4 



where <7j(s) = yT— 4m 2 /s and Qi — 



3 f^o f 2 



[L 3 s - L 5 (m 



TV' 



T+mf/q, 

Ptv+PK° 



2 

max ' 
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FIG. 1. 4> — > tt + iy decay through a photon. 
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FIG. 2. Kaon loop contribution to the p — ui mixing. 
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FIG. 3. Two step mechanism for (j> ^ decay. 
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Fit 


u 


L 2 


L 3 


u 


u 


2L(j + Lg 


U 


Qmax 

(MeV) 


set 1 


0.91 


1.61 


-3.65 


-0.25 


1.07 


0.58 


-0.4 


666 


set 2 


0.91 


1.61 


-3.65 


-0.25 


1.07 


0.58 


0.05 


751 


set 3 


0.88 


1.54 


-3.66 


-0.27 


1.09 


0.68 


0.10 


673 




Li 


L 2 


L 3 


U 


L 5 


21/6 + L$ 


L 7 


A* 


ChPT 
(3l) 


0.4 

±0.3 


1.4 

±0.3 


-3.5 
±1.1 


-0.3 
±0.5 


1.4 
±0.5 


0.5 
±0.3 


-0.4 
±0.2 





TABLE I. 



Different sets of chiral parameters (in 10 -3 units) that yield reasonable fits to the meson-meson 
scattering phase shifts. We have used a hat to differentiate them from those obtained within standard 
ChPT |3l], since in our case we have already differences at the 0(p 4 ) with respect the next-to- leading 
ChPT amplitudes and we have used high energy data in the fit. However, as it is explained in the 
text, we still expect them to be relatively similar once the scales are chosen appropriately (roughly 
fi ~ 1.2q,„ax, see Q for details). 
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FIG. 5. 

Coupled channel IAM results for meson-meson scattering. The dashed, continuous and dotted lines 
are obtained, respectively, with the chiral parameter sets 1, 2 and 3 given in table 1. Note that 
they are indistinguishable for almost every channel. The experimental data for each plot, starting 
from left to right and top to bottom, comes from ||,||, |||^,||, H, 

]§,(§, @,|I§ and, finally [(§ 
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FIG. 7. The (/n 8 , r + 7r - coupling interpreted as a Qs — P mixing and a p — > 7r + 7r decay. 
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